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ABSTRACT 

We observed SAX J1808. 4-3658 (1808), the first accreting millisecond pulsar, in deep quiescence with 
XMM-Newton and (near-simultaneously) Gemini-South. The X-ray spectrum of 1808 is similar to 
that observed in quiescence in 2001 and 2006, describable by an absorbed power-law with photon index 
1.74 ± 0.11 and unabsorbed X-ray luminosity Lx = 7.9 ± 0.7 x 10^^ ergs s-\ for Nh = 1-3 x lO^^ 
cm~^. Fitting all the quiescent XMM-Newton X-ray spectra with a power-law, we constrain any 
thermally emitting neutron star with a hydrogen atmosphere to have a temperature less than 30 eV 
and LAr5(0. 01-10 keV)< 6.2 x 10^° ergs s~^. A thermal plasma model also gives an acceptable fit 
to the continuum. Adding a neutron star component to the plasma model produces less stringent 
constraints on the neutron star; a temperature of 36tg eV and L7vs(0. 01-10 keV)= LS^Qg x 10^^ 
ergs/s. In the framework of the current theory of neutron star heating and cooling, the constraints 
on the thermal luminosity of 1808 and IH 1905+000 require strongly enhanced cooling in the cores of 
these neutron stars. 

We compile data from the literature on the mass transfer rates and quiescent thermal flux of the 
largest possible sample of transient neutron star LMXBs. We identify a thermal component in the 
quiescent spectrum of the accreting millisecond pulsar ICR J00291+5934, which is consistent with 
the standard cooling model. The contrast between the cooling rates of ICR J00291-f 5934 and 1808 
suggests that 1808 may have a significantly larger mass. This can be interpreted as arising from 
differences in the binary evolution history or initial neutron star mass in these otherwise similar 
systems. 

Subject headings: binaries : X-rays — dense matter — stars: pulsars — stars: neutron 



1. INTRODUCTION 

The X-ray transient SAX J1808.4-3658 (hereafter 
1808) has provided many fundamental breakthroughs 
in the study of accreting neutron stars (NSs). It 
was discovered in 1996 by BeppoSAX's Wide Field 
\ Cameras, and type I X-ray bursts were seen, iden- 
tifying it as an accreting NS and constraining the 

■ dista nce ()in 't Zand et al.l1l998t iGallowav fc Cummind 
\ |2006[ ). Coherent millisecond X-ray pulsations, the first 

discovered in accreting systems, were identified during an 

■ outburst using the Rossi X-ray Timing Explorer (RXTE; 
\ IWiinands fc van der Klij|1998[ ). Burst oscillations have 

also been seen at 1808's 401 Hz spin frequency, confirm- 
' ing that thermonuclear burst oscillations in low-mass X- 
ray binaries (LMXBs) represent the spin period o f the 
NS (|in't Zand et al.l[200ll : IChakrabartv et al.ll2003D . 

Transiently accreting NSs in quiescence are seen to 
have a soft, blackbody-like X-ray spectral component, 
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and/or a harder X-ray compon ent generally fit by a 
power-law of photon index 1-2 (jCampana et al.lll998D . 
The harder component is of unknown origin; an ef- 
fect of continued accretio n or a shock from a pulsar 
wind have been suggested (jCampana et al.lll998f ). The 
blackbody-like component is generally understood as 
the radiation of heat from the NS surface. This heat 
is produced by nuclear fusion in the deep crust dur- 
ing accretion, and is radiated from the surface on a 
timescale of lO'* years, prod ucing a steady quiescent ther- 
mal NS luminosity (Brown et al.]ll998l : ICampana et al.l 
ll998HHaensel fc Zdunikil990l ). The deep crustal heating 
rate can be computed if the mass transfer rate is known 
(or estimated). Some transiently accreting NSs have 
been shown to have very low quiescent thermal X-ray lu- 
minosities, indicating either enhanced neutrino emission 
or mass transfer rat es much lower than have been pre- 
viousl y inferred (e.g. IWiinands et ani2001l; Ijonker et al.l 
12004 iTomsick et all 120041 : Honker et al.l l2006f l. The 
coolest of these provide the strongest constraints to 
date on neutrino cooling from NS cores, as a broader 
range of cooling rates is required from X-ray transients 
than from young cqohng pulsars (cf. iPage et ani2004l : 
lYakovlev fc Pethick|[200l . 

X-ray observations in quiescence have shown 1808 
to have one of the lowest quiescent thermal lumi- 
nosities yet mea sured from any accreting neutron star 
(ICampana et all l2002l : iHeinke et all l2007l : iJonker et all 
[2007^. 1808 is of pa rticular import ance because of its 
well-known distance (jGallowav fc Cu mmiug 200(^) and 
relatively stable (within a factor of 2) time-averaged mass 
transfer rate onto the NS measured over multiple out- 
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bursts. 1808's low quiescent thermal luminosity indicates 
that most of the heat absorbed by the NS core during ac- 
cretion is reradiated not as thermal X-ray emission, but 
throu gh neutrino cooling processes (|Yakovlev fc Pethickl 
[2OOI . 

We have obtained a third XMM observation of 
1808 in quiescence in 2007, in conjunction with near- 
simultaneous (separated by 6.5 hours) Gemini optical {g' 
and i') imaging. The science goals included further con- 
straining the thermal component of the X-ray emission, 
constraining X-ray variability in quiescence, and mea- 
suring the sinusoidal orbital optical modulation nearly 
simultaneously with an X-ray observation to determine 
the origin of the optical modulation. One key finding 
is that SAX J1808. 4-3658 provides one of the two most 
constraining upper limits on the thermal component for 
any accreting neutron star. Accordingly, we report on 
the X-ray analysis here, along with comparison to other 
X-ray transients in quiescence, while the optical analy- 
sis and comparison of the X-ray and optical results are 
presented in a companion paper (jPelove et al.ll20d8[ ). 

2. DATA REDUCTION 

We observed 1808 on March 10-11, 2007 (ObsID 
0400230501; starting at 16:24 UT) for 49 ksec with 
XMM's EPIC cani era, using two MOS CCD detectors 
(jTurner et al.ll200ip with med ium filters and one pn CCD 
detector (jStriider et al.ll2001h with a thin filter. All data 
were reduced using FTOOLS and SAS version 7.0.0. Soft 
proton fiares were excluded by excluding times when the 
total MOS count rate exceeded 4 0.2-12 keV counts per 
second, and times when the total pn count rate exceeded 
20 0.2-12 keV counts per second. This left 36.9, 49.7, and 
49.8 ksec in the 2007 pn, MOSl, and M0S2 datasets. 
Event grades higher than 12 were also excluded. We ex- 
tracted spectra from a 10" circle around the position 
of 1808, and combined the pair of simultaneous MOS 
spectra and responses using FTOOLS. We generated re- 
sponse and effective area files using the SAS tasks rmf- 
gen and arfgen, and produced background spectra from 
90" circular source-free regions on the same CCD. The 
spectra were grouped to >15 counts per bin for the MOS 
data, and >30 counts per bin for the pn data (other 
choices gave similar results). 

2.1. X-ray Variability 

We produced background-subtracted lightcurves of 
the 2007 pn data within SAS, and analyzed them us- 
ing HEASARC's XRONOS software ^. Kolmogorov- 
Smirnov and tests on the first 37 ksec of 0.2-12 keV 
pn data (mostly unaffected by background flaring) re- 
vealed mild evidence of variability, as the probability of 
a constant fiux is 3 x 10~^ and 3 x 10"'^ for the two tests 
respectively. 

3. X-RAY SPECTRAL ANALYSIS 

Our X-ray spectral analysis includes photoelectric ab- 
sorption (XSPEC model phabs), with a hydrogen column 
densit y, Nh, fixed at the interst ellar value of 1.3 x 10^^ 
cm-2 (jPickev fc LockmanI fT990l) . We checked this Nh 
value by analyzing a series of Swift observ ations taken 
during the tail of 1808's 2005 outburst (|Konel l2005l : 

^ http: / /hcasarc. gsfc.nasa.gov/docs/xanadu/xronos/xronos. html 



ICampana et al.l [2005af) . We used window timing data 
from June 17, 23, and 29, and July 7 and 13, and photon 
counting data from June 17 and 23. The June 17 obser- 
vation suffered from pilcup in photon counting mode. We 
addressed this by excluding the central 18" of the point- 
spread-function (chosen by fitting a King model to the 
radial profile ^'^). An absorbed power-law model (with 
only Nh fixed between observations) gave a nearly rea- 
sonable fit to the 7 Swift spectra (x^=1.19, null hypothe- 
sis probability=1.4x 10""'), while an absorbed power-law 
plus blackbody model gave a slightly better fit (x^=1.15, 
null hypothesis probability=1.8x 10~^). Photon indices 
ranged from 1.8 to 2.3 for the simple power-law fits, or 
2.2 to 2.4 for the power-law plus blackbody fits, while the 
blackbody temperatures ranged from 0.7 to 1.2 keV. The 
best fit Nh for either fit was 1.2±0.1 x 10^-^ cm~^, consis- 
tent with the Dickey value. We also tested models with 
photoelectric absorption as a free parameter, finding Nh 
consistent with the outburst value. Quoted errors are at 
90% confidence. 

We simultaneously fit our pn and MOS spectra along 
with spectra from the two previous XMM observa- 
tions from March 24, 2001 and Sept. 14, 2006 (see 
iHeinke et al.ll2007i for details of these observations) . The 
X-ray characteristics a re similar to those observed in 
the prior observations (jCampana et al . 2002); the spec- 
trum can be well fit with a power-law of photon index 
1.74 ± 0.11, while a hydrogen-atmosphere model (the 
NSATMOS^^ mod el of (Heinke et all (^OOl . or the sim- 
ilar NSA model of lZavlin etall (|l996( n gives very poor 
fits. 

Although each of the spectra are well-fit by an ab- 
sorbed power-law with the Nh observed in outburst, re- 
quiring all spectral parameters to be identical across ob- 
servations produces somewhat poor fits {xt ~ 1-09, null 
hypothesis probability or nhp=0.03). Allowing a con- 
stant normalization to vary between the observations, we 
find that the 2001 fiux is 1.29t J5 of the 2006 flux, and 

the 2007 flux is lAGt'lj of the 2006 flux. This suggests 
that at its lowest flux levels, 1808 remains X- ray variable 
in quiescence (see also lCampana et aI]l2008a^ ■ Such vari- 
ability could indicate that the X-ray emitting process is 
powered by time-variable accretion, either onto the NS 
or in the interaction of a disk and the NS magnetosphere. 

Allowing only the Nh to vary improves the fit only 
slightly (nhp=0.04), and we note that a different Nh 
in outburst than quiescence is rarely seen in quiescent 
LMXBs, e specially those like 1808 a t relatively low in- 
clination (jJonker fc Nelemand [200I . We fix the Nh 
at 1.3 X 10^^ cm~^ and free the power-law photon in- 
dex and normalizations, including also a NSATMOS 
component with NS mass=1.4 MQand radius 10 km, 
with the same temperature between observations. We 
freeze the distance to 1 808 a t 3.5 kpc, as measured by 
iGallowav fc Cummingl ()2006l ). The parameters of this 
fit are listed in Table 1, and it is shown in Figure [TJ 
No thermal component is required, but a thermal com- 
ponent with kT <30 eV (90% confidence) is permitted, 



Following 'http: / /www. swift. ac.uk/pileup.htmlj 

After this paper was submitted, an indepen dent analysis of 

the Sw ift data found Nh = 1.3±0.1 X 10^^ cm^-^ l|CamDana et al.l 

I2008al) . 

^'^ http: / / xspec.gsfc.nasa.gov/docs / xanadu / xspec / models / nsatmos.html 
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Fig. 1. — XMM X-ray spectra (data and best-fit power-law 
model) of 1808. Solid (blue) line with filled squares: 2007 pn data, 
model. Solid (magenta) line with crosses: 2007 MOS data, model. 
Other data is from the 2001 MOS (black), 2006 MPS (green), and 
2006 pn (red) XMM observations i lHeinke et al.ll2007l) . The max- 
imum acceptable contribution (90% conf.) from the NSATMOS 
component is indicated by a dotted line, prominent at low ener- 
gies. (See the online edition of the Journal for a color version of 
this figure.) 

thus placing a limit on the NS's thermal bolometric (0.01- 
10 keV) luminosity Lnsmi < 6-2 x 10^° erg s~^ This is 
a substantially tighter constrai nt on the NS ' s quie scent 
thermal luminosity than that of lHeinke et al.l ()2007[ ) , and 
may be the tightest constraint on the thermal bolomet- 
ric luminosity of a NS in an X-ray transient. (Note that 
1H1905-I-000 h as a much lower to tal X-ray luminosity 
than 1808; see iJonker et al.l l2007a[) . This tighter con- 
straint is produced by doubling the XMM pn exposure, 
and thus nearly doubling our sensitivity. 

The rather t i ght distance limits of 
iGallowav fc Cummingj (|2006[ ) (3.5 ± 0.1 kpc) produce 
only a 6% uncertainty in our upper limit. Increasing 
the assumed Nh by 10^° cm~^ increases this limit by 
~ 20%; decreasing the Nh decreases the limit similarly. 
Changing the assumed NS mass to 2 M0(as suggested 
by the results of Deloye et al. 2008) increases the upper 
limit by 10%, as does altering the assumed NS radius to 
13 km. 

The nature of the high-energy spectral component in 
qLMXBs is unclear, and thus other continuum mod- 
els are possible. One possibility for the emission is a 
hot plasma, possibly a shock between an infalling ac- 
cretion stream and a pulsa r wind (e.g. ICampana et al.l 
[2002: Bogd anov et al.ll2005l ). Such a model is not ruled 
out by the spectral fits, which allow (x^ = 0.95 for 84 
degrees of freedom) a hot plasma MEKAL model with 
kT =4-8 keV instead of a power-law spectrum. Adding 
a faint NSATMOS component is modestly preferred in 
this spectral fit; the best fit has kT = 36lg eV, for a 

bolometric NS luminosity of Lx = l-3to8 ^ ^^^^ crgs/s 
(Table [T]). This model allows a larger NS component due 
to the flatter shape of the model at low energies, and 
demonstrates the dependence of our constraints on the 
continuum model. The nature of the hard spectral com- 
ponent, and thus choice of MEKAL vs. power-law model, 
could be tested by a sensitive search for line emission, or 
by measuring the hard tail out to higher X-ray energies, 
with next-generation X-ray instruments. 



4. COMPARISON TO OTHER RESULTS 

iHeinke et al.l (|2007t ) noted a possible trend of NSs with 
low mass transfer rates having faster cooling. To test 
this, we have increased our sample of LMXBs with use- 
ful constraints on their mass transfer rates and quiescent 
thermal luminosities (Table [2]). For several systems, we 
take both values from the recent literature. For some, 
we compute the quiescent 0.01-10 keV thermal NS lumi- 
nosity given the 0.5-10 keV fluxes in the literature and 
NS atmosphere models, or by reanalyzing Chandra or 
XMM data on quiescent NSs. For several systems, we 
estimate the time-averaged mass transfer rate (or an up- 
per limit) from X-ray flux histories in the literature, since 
the RXTE ASM lightcurvcs (used in iHeinke et anl2007t l 
generally do not cover their full outbursts. Distance es- 
timates are often rather poorly quantified. We take best 
estimates of these distances from the quoted references, 
noting that an uncertainty of 50% in distance is a reason- 
able upper bound. Changes in distance at this level will 
not greatly affect our (broad) conclusions. We also list 
orbital periods where known. We discuss a few systems 
in detail below. 



4.1. Comparison to other LMXBs in quiescence 

iJonker et al.l ()2007a[ ) obtained a long observation of the 
quiescent NS LMXB lH1905-f000, deriving a very tight 
limit on the quiescent bolometric NS flux of < 1 x 10^^ 
crgs/s (for a 0.1 keV blackbody). The mass transfer rate 
of 1H1905-I-000 is not known; but for any mass transfer 
rate > 10^^^ Mq/jx, the quiescent l uminosity limit is be- 
low th e fiducial pion cooling curve of lYakovlev fc Pethickl 
(|2004[ ). On the other hand, the mass transfer rate for 
1808 is rather well-known, as the mass transfer rate can 
be inferred from RXTE ASM monitoring of multiple out- 
bursts, and by predictions fro m gravitational radiation 
(iBildsten fc Chakrabartvl l200lD. which agree very nicely. 
IJonker et al.l ()2007bD used Chandra to identify the 



likely quiescent X-ray counterpart to 1M1716-315, with a 
soft thermal spectrum. W e assume H-poor material, as 
suggested bv IJonker et a l.' (2007b!), to select the larger 



distance of IJonker et al. (.2004)) . We ut ilize the histori- 



cal X-ray fluxes of I Jonker "eTall (I2007br i to estimate an 

average mass transfer rate of M~ 2.5 x Mo/year 
o ver 37 years, wh i ch we take as an upper limit. 

iTomsick et al.l (120071 ) identified the quiescent X-ray 
co unterpart to 4U17 30-22. We use the outburst history 
in lChen et al.l (|l997l ) to compute its time-averaged mass 
transfer rate over the past 37 years. 

4U 2129-1-47 was a bright X-ray source from 1971 until 
sometime between 1980 and 1983, and has been in qui- 
escence since then (jPietsch et al.lll986D . 4U 2129-1-47 is 
optically bright enough {B ~ 16.8 — 18.3) that historic 
photographic plates are useful; they show 4U 2129-f47 
to be "on" from 19 63 to 1979, and "off" from 1938 to 
1943 (|Wenzel|[l98l . We estimate an upper hmit for 
the average X-ray flux of Fx ~ 8 x 10~^° ergs/s (0.1- 
20 keV), maintained over 20 of the last 45 years. Us - 
ing a distance of 6.3 kpc (jCowlev fc Schmidtkd Il990l ). 
we derive an upper limit of Lx < 4 x 10"^^ ergs/s. How- 
ever, the lightcurve in outburst shows extended V-shaped 
eclipses, indicating that we observe only scattered light 
from an accretion disk corona. This may decrease our ob- 
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served flux, possibly by a factor of more than 10 P The 
peak lumino sity of a non-radiu s -expan sion X-ray burst 
observed by iGarcia fc Grindlavl (|1987[ ) is i = 5 x lO^^ 
ergs/s. This is a factor of 40 lower than the Eddington 
limit for hydrogen-rich material. Allowing for a factor of 
40 reduction in our observed flux, we estimate a (con- 
servative) upper limit of Lx < 1-6 x 10^® ergs/s, and 
M< 5.2 X 10"^ Mo/year. 

4.2. Comparison to other accreting MSPs in quiescence 

Of the other accreting millisecond X-ray pulsars, deep 
quiescent studies have not been performed for Swift 
J1756.9-2508 and HETE J1900. 1-2455. For most of 
the remainder, we compute mass transfer rates using 
the distance estimates, time-average d bolometric flu - 
ences, and recurrence time limits in IGallowavl (|2006f ). 
For XTE J1 814-338, we assume a rec urrence time 
of 19 years (jWiinands fc Revnoldsl I2003D . For XTE 
J 175 1-305, the discovery outburst seems to have been 
much brighter than the other three recorde d outbursts 
(|Markwardt et al.ll2002i : iLinares et 311120071 ). so we add 
the estimated outburst fluences and average over the past 
12 years. 

XTE J1814-338 (1814) docs not have a pubhshcd qui- 
escent flux measurement, but one archival (Sept. 6, 2005) 
XMM observation exists. We have analyzed this dataset 
to search for 1814 in quiescence. We excise periods of 
high background, >30 (> 6) cts/s in the 0.2-10 keV pn 
(MOS) data, giving 22.9, 30.9, 29.5 kscc of pn, MOSl, 
and M0S2 data respectively. A slight flux enhancement 
may be seen in the 0.2-4 keV images, although it is not 
highly significant. We estimate the flux of 1814 by com- 
puting the counts within 15" of 1814's nominal position, 
and subtracting background estimated from an annulus 
from 20" to 40" away. Calculating the fluxes for an ab- 
sorbed power-law of photon index 2 using PIMMS gives 
an (averaged) (0.5-10) of 6.3 ±4.3 x IQ-^^ ergs cm'^ 
s~^. Using a blackbody of temperature 0.2 keV gives 
Fx = 3.9 ± 1.6 X 10~^^ ergs cm^ s~^. We therefore iden- 
tify Fx < 9 X 10~^^ ergs cm~^ s~^ as a 3cr upper limit on 
the blackbody-like flux^^ from a NS in 1814, and for an 
assumed distance of 8 kpc, LArs(0. 01-10 keV)< 1.7x 10^^ 
ergs/s. 

IGR J00291-h5734's (hereafter 00291) quiescent X-ray 
flux has bee n observed wi t h Cha n dra on five oc c asions , 
repo rted bv iJonker et all (|2005D . iTorres et all (|2008D . 
and iJonker et al.l (|2008f l. The third observation found 
00291 to be significantly brighter (factor of two), with 
a softer spectrum, than th e flrst two or the later obser- 
vations (jJonker et al.ll2005l ). This observation can be fit 
with a NS atmosphere model alone, but the implied ra- 
dius is rather small (only 3'^i km for a distance estimate 
of 4 kpc; see Galloway et al. 2005, Torres et al. 2008, 
Jonkcr ct al. 2005). However, the addition of a power- 
law component, dominant at high energies, allows a rea- 
sonable NS radius. The temperature of the NS appears 
higher in this observation than in other observations. Ex- 
ploring the nature of this variation lies outside the scope 




IJonker fc van der KlisI 1)20011 ') estimate that the luminosity of 
the accretion disk corona source 2A 1822-371 is reduced by a factor 
of 40 due to its high incUnation. 

^'^ Since a harder spectrum gives a higher upper limit, the total 
luminosity could be slightly higher. 
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Fig. 2.— X-ray spectrum of IGR J00291-I-5934, fitted with a 
power-law alone. Data and model in upper panel, residuals in 
lower panel. 

of this work. Since only the faintest blackbody-like flux 
measurement can represent the emission from a slowly 
cooling NS heated by multiple outbursts, this observa- 
tion is in any case not relevant to our purpose. 

The other four observations of 00291 show a similar 
spectrum and X-ray flux (of 7 x lO"^** ergs cm~^ s~^, 
0.5-10 keV, for a power-law fit). We analyze the longest 
observation, from 2005 November (it has 143 counts, vs. 
36, 63, and 22 from the othe rs), with shg h tly di fferent 
methods than those used by iTorres et al.l (|2008( ) . We 
fit a model consisting of an absorbed powerlaw, with an 
optional NS atmosphere, with Nh fixed to 4.6 x 10^^ 
(the average value from several methods summarized by 
ITorres et al.ll2008l ). We binned the source spectrum so 
that each bin contains at least 10 counts (using 15 or 
20 counts/bin y ields similar results) . We set the NS dis- 
tance to 4 kpc (jGallowav et al.l [20051 ). radius to 10 km, 
and mass to 1.4 Mq. We find that the power-law compo- 
nent is absolutely required, but that the NS component 
significantly improves the fit (see Table |3l Figures |2| and 
|3|). The fraction of the 0.5-10 keV flux attributable to 
the NS component is 44^11%- 

The reduced is reduced by half with the addition 
of the NS component, and the temperature of the NS 
component is inconsistent with the model minimum (ef- 
fectively zero) at the >99% level (Ax^=7.0). An F-test 
indicates that the probability of attaining an equivalent 
fit im provement by adding on e model parameter is only 
0.4%. iProtassov et all (|2002D showed that the F-test is 
often inaccurate for testing the necessity of adding an ad- 
ditional spectral component. We therefore simulated 100 
datasets using the best-fit absorbed power-law model, 
and fit them with this model and with the absorbed NS 
plus power-law model. None of our simulations gave 
a Ax^ or F-statistic larger than that produced by our 
model, allowing us to conclude that the probability of 
concluding incorrectly that a NS component is required 
for 00291 is less than 1%. 

The difference between our results and those of 
ITorres et all ()2008[) (which found an upper limit of 19% 
to a thermal component for 00291) may be attributed to 
four factors: our inclusion of the lowest energy bin in our 
fit (this bin is dominated by data above 0.5 keV), our use 
of a fixed distance (which imposes a relation between NS 
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Fig. 3.— X-ray spectrum of IGR J00291+5934, fitted with a 
power-law and NS atmosphere model. The NS atmosphere com- 
ponent is shown by the dashed line, and the power-law component 
by the dotted line, in the upper panel, while residuals are shown 
in the lower panel. Note the improvement in the fit compared to 
Figure E] 

temperature and NS flux), Torres et al's use of a fixed 
NS temperature, and, perhaps most importantly, Torres 
et al.'s method of first fitting an absorbed power-law, 
and then fixing the best-fit power-law parameters before 
adding and constraining a blackbody component. This 
method is likely to underestimate the flux that may be 
present in a second component, since the parameters of 
the first component are not allowed to vary from those 
of the best one-component fit. We have replicated the 
spectral fit of Torres et al., with similar results (the NS is 
< 19% of the 0.5-10 keV unabsorbed flux), but when we 
free the power law component parameters the NS com- 
ponent then makes up 56^35% of that flux. Torres et al's 
method is co mmon in the literature fe.g. lCampana et al.l 
[200ll2005br i. but we do not feel it is the most appropri- 
ate when constraining the contributions of broad band 
components. 

The best estimates of the mass transfer rates and qui- 
escent X-ray luminosities in Table[2]are plotted in Figure 
[31 These estimates suffer uncertainties in distances and 
recurrence times (where appropriate, we plot upper lim- 
its, in some cases upper limits in both mass transfer rate 
and Lj^s)- However, there is clearly a large variation 
in cooling behavior exhibited by these accreting neutron 
stars, with quiescent therr nal luminositi e s rang ing by sev- 
eral orders of magnitude. iHeinke et al.l (|2007| ) suggested 
that NSs with lower mass transfer rates, which are gen- 
erally older systems, have experienced significant mass 
transfer over their lives and are now more massive, and 
inclined to faster cooling. The larger dataset collected 
here argues against this idea. In addition to the two 
systems with high mass transfer rates, 4U 1730-22 and 
00291 lie near the predictions of standard cooling. We 
note the contrast between the bright thermal component 
of 00291 and the faint thermal component of 1808. These 
two systems have similar periods (2.46 and 2.01 hours, 
respectiv ely) and heated low-mass (brown dwarf) com- 
panio ns (jBildsten fc Chakrabartvl[200ll : iGallowav et al.l 
l2005f ). One might expect similar evolutionary histories, 
similar amounts of accreted mass, and thus similar cool- 
ing rates from these two systems. The observed dif- 
ference in cooling requires different neutrino emission 
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Fig. 4. — Measurements of, or limits on, the quiescent ther- 
mal luminosity of various NS transients, compared to estimates of, 
or upper limits on, their time-averaged mass transfer rates. Data 
from the compilations of IHeinke et al.l 1)20071 1 and Table 3. The 
predictions of standard NS cooli ng and various enhanced c ooling 
mechanisms are plotted, foUowing lYakovlev fc Pethickl l l200l'l . The 
accreting millisecond pulsars (including the intermittent ones) are 
labeled separately (in red). The effect of using a MEKAL rather 
than a power-law continuum for SAX J1808. 4-3658, allowing a 
higher upper limit on the NS component, is indicated with the 
dotted line. The effect of increasing the distance by a factor of 1.5 
for any system is indicated with an arrow (labeled Dxl.5). 

mechanisms, suggesting that 1808's NS may be signif- 
icantly more massive than that in IGR J00291-I-5934. 
This might be interpreted as due to differences in their 
mass transfer histories, or by the NS in 1808 being born 
with a higher mass. 

There are a few caveats to consider. The evidence for 
the thermal component in 00291 is only at the 3a level. 
It is not absolutely certain that the thermal component 
in 00291 is produced from deep crustal heating, espe- 
cia lly considering the o dd quiescent X-ray behavior seen 
bv lJonker et al.l (f2005[ ). The distance and mass transfer 
rate of 00291 are not known to high accuracy. However, 
changes in 00291 's distance by a factor of two, or mass 
transfer rate by a factor of 10 or more, would not elimi- 
nate the contrast between 1808 and 00291 (see Fig. 4). 

5. CONCLUSIONS 

The three combined XMM observations of SAX 
J1808. 4-3658 in quiescence provide one of the most strin- 
gent constraints on the thermal component of any tran- 
siently accreting neutron stars observed so far, with the 
thermal L^sMi < 6 x 10"^° ergs/s (for a powerlaw con- 
tinuum model) or < 1.9 x 10^^ ergs/s (for a MEKAL 
thermal plasma model). Combined with 1808's well- 
constrained mass transfer rate and distance, this con- 
straint strongly requires enhanced neutrino cooling from 
the NS in 1808. The models of kaon and pion cool- 
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ing presented bv lYakovlev fc Pethickl ()2004D are excluded 
by this constraint, favoring direct URCA neutrino emis- 
sion processes involving protons, hyperons or deconfined 
quarks, depending on the constituents of matter at supra- 
nuclear densities. 

We have compiled literature estimates of the mass 
transfer rates and quiescent thermal luminosities of a 
large number of transient NS LMXBs. Uncertainties in 
transient outburst histories and the lack of detections 
of thermal components limit the number of useful data 
points. Although many of the measurements are upper 
limits, there is strong evidence for a range of NS neutrino 
cooling rates. Uncertainties in the distances, while sig- 
nificant for many systems, affect quiescent and outburst 
flux equally (see Fig. 4), and thus do not greatly affect 
these results. 

The accreting miUisecond pulsar IGR J00291-I-5934 
shows evidence for a thermal component in its quies- 
cent X-ray spectrum, which is consistent with the pre- 
diction from standard cooling. 00291 seems similar to 
1808 (both show heated brown dwarf companions, and 
expanding ^2-hour orbits). Thus 0029 1's much hotter 
NS (implying a less massive NS) comes as a surprise. 
If 0029 1's thermal emission is a result of deep crustal 
heating, then it is likely that 00291 either experienced a 



different binary evolution than 1808 (e.g. starting with 
a smaller-mass companion?), or began its evolution with 
a smaller- mass NS. 

Note added in proof- ICampana et al.l (|2008bl submit- 
ted) have independently found evidence for a thermal 
component in 00291 from XMM data, confirming our 
detection of such a component. 
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TABLE 1 

Spectral Fits to SAX J1808.4-3658 



Epoch 


Nh X 10^1 


F/kT 




NS kT, eV 


LNS,bol 




Power-law -|- NSATMOS Fit 


2001 
2006 
2007 


(1.3) 


1 o,+0.16 
1-»3_0.17 
1 74+0.11 
1- '^-0.12 


7.6+14 X 10^^ 0.83/83 
5.lt° j X 10=^1 
7.9±0.7 X 10^1 


< 30 


< 6.2 X 10^0 


MEKAL + NSATMOS Fit 


2001 
2006 
2007 


(1.3) 


0+21 
5 4+1-9 


7.5+^-2 X 10^1 0.90/83 
5.110;? X 10=^1 
7.6±0.7 X 10^^ 


361^ 


1.31";^ X 10^1 



Note. — Simultaneous spectral fits to SAX J1808. 4-3658 data, from 2001, 2006 and 2007, to either a power-law + 
NSATMOS model, or MEKAL + NSATMOS model. Parameters and results from the fit as a whole are listed only in 
the first line for each fit. Errors are 90% confidence for a single parameter. A'^^ is held fixed at the interstellar value, 
in agreement with Swift measurements in outburst (see text). Luminosities in erg s~i; Lx for 0.5-10 keV, Lj^gi^^i 
for 0.01-10 keV. 



TABLE 2 

Estimates of Quiescent NS Luminosities and Mass Transfer Rates 



Source 


^orh 


Nh 


kT 


D 


Outbursts 


Timeline 


M 


Lns 


Refs 




(Hours) 


(10^1 cm-2) 


(eV) 


(kpc) 




(years) 


(A%yr-^) 


(erg s-i) 




SAX J1808-3658 


2.01 


1.3 


< 30 


3.5 


5 


12 


9 X 10"^2 


< 4.9 X 10^0 


1,2,3 


IH 1905+000 


< 1.5? 


2.1 


<40 


10 


1 


34 


< 1.1 X lO^io 


< 1.0 X 10^1 


4,5 


2S 1803-45 


~ 9? 


1.47 


<91 


7.3 


2 


33 


< 7 X 10-" 


< 5.2 X 10^2 


6,2 


1M1716-315 


< 1.5? 


7.4 


116 


6.9 


1 


37 


< 2.5 X lO^io 


1.3 X 10^3 


7,1 


4U 1730-22 




3.7 


131 


10 


1 


37 


< 4.8 X 10-" 


2.2 X 10^3 


8,9,1 


4U 2129+47 


5.24 


1.7 


119 


6.3 


1 


45 


< 5.2 X 10-9 


1.5 X 10^3 


10,11,12,1 


XTE J1751-305 


0.71 


9.8 


< 71 


8 


4 


12 


6 X 10-^2 


< 4 X 10^2 


13,14,15,1,2 


XTE J0929-314 


0.73 


0.76 


< 50 


10/9 


1 


12 


< 2.0 X 10-" 


< 1.0 X 10^2 


16,17,18,15,19 


XTE J1807-294 


0.67 


4.6 


< 51 


8 


1 


12 


< 8 X 10-12 


< 1.3 X 10^2 


20,17,19,2 


XTE J1814-338 


4.27 


1.6 


< 69 


8 


2 


19 


3 X 10-12 


< 1.7 X 10^2 


21,22,17,2 


ICR 00291+5934 


2.46 


2.8 


71 


4 


3 


12 


2.5 X 10-12 


1.9 X 10='2 


23,24,25,2 



Note. — Estimates of quiescent thermal NS luminosities from X-ray transients, orbital periods (those followed by ? are rough estimates- 
see the references), assumed distances, and mass transfer rates (inferred from RXTE ASM observations for systems with RXTE-era outbursts, 
best historical data otherwise). Quiescent thermal luminosities are computed for the unabsorbed NS component in the 0.01-10 keV range. 
Outbursts and years columns give the number of outbursts seen and the time baseline used to co mpute M. References as follow s: 1 : Mass trans- 
fer rat e co mputed in this work, 2: Q uiescent NS thermal lum inos ity computed i n this work . 3: IGalloway &: C ummind ll200i). 4: IJonker et all 
1I2OO6I). 5: IJonker et all Il2007all. 6: IC ornelissejTan (12003 ), 7: Donker et all || 2p07b( ), 8: ITomsick e t alJ l|2007l), 9: IChcn et alj l|1997|) ,l^ 
INowak et aLI II2002I). 11: IPietsch et al.l (1986). 12: Wenzer il983l . 13: Markwardt et al. (2002), 14: Miller et al. (2003), 15: Wiinands c TaLl 
(120051). 16: IGalloway et al.l (120021). 17: iGallowav, (,2006,). 1 8: IJuetTeral., (.200a) . 19: iCampana et al.. (,2005b.). 20: .Canipana fc Stella. (.20031 ). 
21: IKrauss et al.l (I2005i )~22: IWiinands fc Revnoldsl (120031) . 23: IGallowav et al J ( 120051) . 24: IJonker et al.ll2005l) . 25: ITorres et al.l (120081) 



TABLE 3 

Spectral Fits to IGR J00291+5934 



Model 


Nh X 102i 


F 


Lx 




kT, eV 


Lnsm 


PL 
NS+PL 


(4.6) 
(4.6) 


2 6+0-5 
1.5+0.7 


1.3 X 10^2 
1.8 X 10=*2 


1.08/12 
0.54/11 


'1-9 


1.9l0;« X 1032 



Note. — Spectral fits to 2005 Chand ra data on IGR J00 291+5934. A distance of 4 kpc is assumed, and 
Nh is fixed to the value determined bv ITorres et al.l ( 120081 ). Luminosities in erg s— i; Lx for 0.5-10 keV, 
LNS.bol f°r 0.01-10 keV. 



